Submarine Landslides and Local Tsunami Waves (Corinthos Gulf, Greece) 



L.I.Lobkovsky^, G.A.Papadopoulos^ ^ I.A.Garagash^^ O.R.Kozyrev^^ and R.Kh.Mazova^ 
(a) Institute of Oceanology, Russian Academy of Science, Moscow, Russia 
(b) Institute of Geodynamics, National Observatory of Athens, Athens, Greece 
{c) Institute of Physics of Earth, Russian Academy of Science, Moscow, Russia 
(d) Nizhny Novgorod State Technical University, Nizhny Novgorod, Russia 

It is considered a problem of submarine sediment slide which generates the surface water waves. To 
simulate numerically the landslide motion it is used the method which permits to take into account 
detailed rheological properties of slide body constituents. The numerical simulation of landslide- 
generated surface water waves is performed on the basis of nonlinear shallow water equations. It 
is studied a landslide behavior for two values of maximal friction angle. The detailed comparison 
of landslide dynamics and evolution of landslide-generated surface water waves during sliding is 
performed. Also, the evolution of dipolar water wave generated in the beginning of sediment sliding 
is studied. It is obtained that this dipolar wave is then transformed in two wave groups: crest and 
trough coming seaward and trough accompanied by crest coming to the shoreline. The second wave 
group leads firstly to sea recession from the beach and only then to large runup (tsunami wave with 
first negative phase). 



Large water waves produced by submarine land- 
slides were observed in many regions of the world 
for last 50 years (local tsunamis) ( see, e.g. [0). 
For numerical simulation of landslide-generated 
tsunami there were proposed a number of mod- 
els, main of which are a rigid-body model and 
viscous (visco-plastic) model (for review, see, e.g. 
Q). First model, because of its specifics, over- 
estimates the water surface response to the sub- 
marine perturbation while second model underes- 
timates it. To simulate adequately the landslide- 
induced tsunami it is necessary to use methodics 
which take into account the both detailed struc- 
ture of landslide body and mechanical characteris- 
tics of slide-body constituents during sliding (see, 
e.g.[^). The character of disconsolidation of con- 
stituents in sediment surface layer appears to be 
a key factor which controls the localization pro- 
cess of the plastic strain and thus the slope insta- 
bility. The model numerical simulation of sliding 
process at the continental slope was carried out for 
a number of slope parameters. In particular, such 
simulation of landslide induced tsunami was per- 
formed for parameters of Corinthos Gulf, Greece 
(see fig.l, 2) where in 7 February 1963 it was ob- 
served a damaging tsunami wave formed without 
any seismic events: mass of unconsolidated sedi- 
ments slumped into the sea water around of the 
local river [Q. It is considered layer-by-layer slid- 
ing of upper part of elastically-plastic sediment 
layer on the slope surface which is formed dur- 
ing the landslide process. It is proposed a dis- 
tinct interface between water and landslide body 
with water density pw{x^z) = const ^ and landslide 
density is a function of coordinate ps = ps{x^z) 
(see, fig.l). The coordinate origin is taken at the 
coastline, with x-axis at the undisturbed seawa- 
ter level and directed seaward and with z-axis di- 
rected upward. In the fig.l z = —h{x^t) is a vari- 




FIG. 1: Model scheme for simulation of tsunami wave 
generated by moving submarine landslide 



able depth of seawater, z = —h{x^O) = —hs{x^O) 
is a submarine slope profile before initial action 
{t = 0). D{x^t) is a landslide body thickness: 
h{x^t) = hs{x^t) — Dix^t). The shape of seawa- 
ter surface z = r]{x^t), ^(^,0) = 0. At simu- 
lation it was used explicit finite-difference scheme 
which permits to simulate the nonlinear behavior 
of pore saturated sediments under conditions of 
plastic flow above yield stress. Also, possible ini- 
tial action at the process (adequate to some seis- 
mic action) was taken into account (see fig. 2, where 
it is presented a geometrical scheme for numerical 
simulation of landslide (fig. 2a) under the action of 
earthquake with moderate size magnitude (equiv- 
alent to slumping volume 57.000m^ in Corinthos 
Gulf), fig.2b). 

It was used a layered model of sediments rested 
on relatively rigid base. For each layer there were 
taken a layer density, shear modulus, bulk mod- 
ulus, cohesion, maximal friction angle, and ten- 
sile strength. At the first stage, it was made the 
simulation of initial, preliminary stressed state of 
slope with sediments formed under action of its 
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FIG. 2: a) Simulation scheme of submarine slope; up- 
per layer is an elastically-plastic sediment mass rested 
on elastic base; b) Change of base velocity during 
earthquake with duration 6sec and frequency 3Hz; 
maximal velocity is equal to 0.2m/ sec corresponding 
to earthquake magnitude equal 6 
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own weight and saturation with the water under 
sea pressure. At the second stage, it was taken into 
account the rheological effect of decrease sediment 
mass strength above yield stress, pore saturation 
of sediment mass and possibility of landslide-mass 
liquefaction under conditions of a seismic or any 
external action. And finally, it was performed nu- 
merical simulation of surface water waves (evolu- 
tion and runup) on the basis of nonlinear system of 
shallow water equations with using of constructed 
explicit difference scheme with fulfilled stability 
conditions. Dynamical interaction between land- 
slide motion and surface waves was taken into ac- 
count via continuity equation. The results are con- 
sistent with those in conventional (rigid-body and 
viscous-fluid) model and agree qualitatively with 
natural data. So, it was really observed the sea 
recession (first negative tsunami wave) in the re- 
gion of estuary of the local river in Corinthos Gulf, 
Greece where coastal strip have been slumped into 
the sea water (see above). 

In the work there are presented results of numer- 
ical simulation for two values of maximal friction 
angle (j)s • 20^ and 32^ (see below). In fig. 3a, b 
there are presented a picture of successive posi- 
tion of landslide surface for several times with 10 
sec interval (dashed curves) calculated for maxi- 
mal friction angle (j)s- 20^ and 32^, respectively. 
Solid line corresponds to initial landslide shape 
{t = 0). At given values of parameters in the 
case (j)s = 20^, landslide comes only finite dis- 
tance, equal to 1200m. The duration of sliding 
process may be determined from the time instant 



FIG. 3: Evolution of submarine slope surface shape 
during landslide movement: a) = 20°; b) (j)s = 32°. 
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FIG. 4: Shear strain distribution at (j)s = 20° a) t = 
10 sec; b) t = 40 sec. 



when velocity of forefront reaches zero, and is near 
50 sec. Maximal velocity of landslide forefront in 
this case is reached at near 40 sec and its value is 
near 6m /sec. In contrast, rear part of the landslide 
shifts to lower x. In the case (j)s = 32^, landslide 
comes only near 750 m (fig. 3b) and its maximal 
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velocity is essentially lower though it is reached 
also at near 40 sec. 

In fig. 4a, b there are presented a shear strain 
distribution in landslide body for maximal friction 
angles (j)s = 20^ for time 10 sec and 40 sec, re- 
spectively. From numerical results, it's seen that 
in the case <ps = 20^ landslide forefront reaches 
level of horizontal bottom, and horizontal region 
near shoreline disappears so that landslide is plane 
up to shoreline. In the case of (f)s = 32^ land- 
slide don't reach a horizontal bottom level so that 
landslide shape is plane only in the middle part of 
slope, and horizontal region near shoreline persists 
(see, 0). 

For numerical simulation of landslide-generated 
surface water waves it was used a system of shallow 
water nonlinear equations with taking into account 
the friction and Coriolis force 
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where x, z are real space coordinates, t is the time, 
u{x^t)^v{x^t) are velocity vector components, 77 is 
the surface elevation amplitude relative to undis- 
turbed state, ^ is a gravity acceleration, H{x^ t) = 
h{x^ t) -\-r]{x^ t) is a total depth, (h is a depth mea- 
sured from undisturbed seawater level), r is a bed 
friction coefficient, / is a Coriolis force parameter. 

Surface water waves were in fact generated by 
moving submarine landslide because of continuity 
equation in the form 
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To solve the system numerically it was used a 
scheme constructed in analogy with those used in 
one-dimensional case. Numerical simulation pro- 
cess was based on splitting of difference operator: 
equations on x and z are integrated separately at 
two semi-steps in time. It was used a space dis- 
persed pattern for different presentation of vari- 
ables. Th avoid the appearance of numerical in- 
stability problem and necessity to use a filtration 
scheme it was applied a method with first order in 
time scheme. The Cartesian coordinate system is 
taken so that total depth is positive in the region 
of liquid phase and is equal to zero at the phase 
interface {r] = —h). So, negative values of H indi- 
cate to solid phase (beach) . To localize the bound- 
ary it was used a linear extrapolation procedure. 
Since it is known that upstream oriented schemes 
are usually stable then it was used a conventional 
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FIG. 5: Surface water wave generation by submarine 
landslide (maximal friction angle is 0s = 20°). 



procedure for advective terms of equations of mo- 
tion and continuity. The procedure is based on 
choosing of direction of space discretization. 

In fig. 5 it is presented an evolution of surface wa- 
ter waves (upper curves) generated by moving sub- 
marine landslide (lower curves) for the case of max- 
imal friction angle = 20^. The numbers near 
upper curves (surface waves) correspond to those 
near lower curves (position of landslide). There 
arise two distinct wave groups: crest and trough 
moving seaward, and second trough moving to a 
beach. It is seen that wave crest velocity is notice- 
ably higher than that of landslide forefront. The 
picture is consistent with results of numerical sim- 
ulation in conventional models (rigid body, viscous 
fiuid, visco-plastic fiuid). 

In fig. 6 there is presented an evolution of 
landslide-induced surface water waves correspond- 
ing to maximal friction angle = 20^ with time 
step 10 seconds. In the beginning, it is generated 
a dipolar wave with trough oriented to a beach. 
This dipolar wave moves in the deep water derec- 
tion. But soon an additional crest in the region of 
the trough appears which moves to a beach. Such 
recession of the sea water indeed was observed in 
the region where landslide occurs which phenom- 
ena leads to anomalous tsunami with first negative 
phase [||, (depressive wave). 

In conclusion, in contrast to kinematic method, 
complete solution of problem with stress strained 
state of slope at seismic action gives the possibil- 
ity to estimate the distribution of residual strain 
and displacement in overall sediment volume with 
taking into account its possible liquefaction and 
disconsolidation Numerical code FLAC, in 

contrast to finite element method realizes an ex- 
plicit finite difference scheme for solution of three- 
dimensional (3D) problems in continuum mechan- 
ics that permits one to simulate nonlinear behavior 
of pore saturated mass under conditions of plastic 
ffow above yield stress. For calculations, material 
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FIG. 6: Evolution of surface water waves generated by 
submarine landslide (maximal friction angle is (f)s — 
20"). 



was divided to polyhedral elements in the limit of 
grid corresponding to shape of calculated object. 
Each element behaves itself according to action of 
applied forces and edge's restrictions. The grid is 
frozen in material and moves together with it be- 
ing undergone to finite strain and displacement. 
Explicit Lagrange scheme of calculations guaran- 
tees an accurate simulation and material flow. A 
key factor to control the process of localization 
of plastic strain and related instability is a char- 
acter of constituents composing the surface sedi- 
ment layer. The essence of this phenomena is that 
straining of real sediment layer after reaching of 
destruction point is continued at decreasing stress, 
i.e. it occurs decrease of limit of strength at in- 
crease of strain up to stress reaches some finite 
or residual level. Process of sediment mass slid- 
ing is essentially determined by friction angle. Ac- 
cording to experiments, submarine sediments are 
disconsolidated when reaching maximal strength. 
Firstly, sediment mass is strained elastically but 
when limit condition r = = ccos(ps — asin(j)s is 
fulfilled it begins to move. In result, friction angle 
decreases from the peak value ^ = (^g in the maxi- 
mal strain point, and it is established a new equi- 
librium state at the level r = Tk = ccos(j)k — o'sin(j)k 
at residual value (j) = (j)k §. Moreover, strength 
decrease of sediment mass during the development 
of plastic strain at static action is a key factor to 
control a slope stability. In the case of neglecting of 
the effect of sediment mass strength decrease slope 
displacements remains to be small which fact leads 
to mistake conclusions. And for further agree- 
ment of numerical simulation results with natural 
data it is necessary to attract the detailed data 
about inner structure of sediment formations in 
landslide dangerous regions of slopes and landslide 
constituent properties. 
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